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1. Introduction 

Recently there has been considerable interest in the novel suggestion that we live 
in a Universe that possesses more than four dimensions. The standard model fields 
are assumed to be confined to a hyper-surface (or 3-brane) embedded in this higher 
dimensional space, in contrast to the gravitational fields which propagate through 
the whole of spacetime [1-10]. In order for this to be a phenomenologically relevant 
model of our universe, standard four-dimensional gravity must be recovered on our 
brane. There are various ways to do this, the most obvious being to assume that the 
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extra dimensions transverse to our brane are compact. In this case gravity can be 
recovered on scales larger than the size of the extra dimensions [5-7]. This is different 
from earlier proposals since the restrictions on the size of the extra dimensions from 
particle physics experiments no longer apply, as the standard model fields are confined 
to the brane. The extra dimensions only have to be smaller than the scale on which 
gravity experiments have probed, of order 0.1mm at the time of writing. Another 
way to recover four-dimensional gravity at large distances is to embed a positive 
tension 3-brane into an AdSr, bulk [9,10]. In this scenario four-dimensional gravity 
is obtained at scales larger than the AdS radius. Randall and Sundrum showed that 
this could produce sensible gravity even if the extra dimension was not compact. 

Several aspects of the above extra dimensional scenarios have since been inves- 
tigated, and compared with the standard four-dimensional case. The cosmology of 
a 3-brane in a five-dimensional bulk was studied and its Friedman equation derived 
and shown to contain several extra terms [11-15]. Pertubations to this homogenious 
case have been examined [16-19] as have some inflationary models [20,21], as well 
as phase transitions, topological defects and baryogenesis [22]. More recently it was 
shown how to embed the Randall-Sundrum models within supergravities [23-25] and 
then within string theory compactifications as in [26-29]. The seemingly arbitrary 
feature of having an AdS bulk spacetime is actually well motivated as it is found 
as a supersymmetric vacuum to supergravity theories, inspiring several more recent 
brane world models [30-36]. Generalisations have been made to six dimensions, and 
it has been demonstrated that both singular and non-singular self-gravitating string 
defects can localise gravity [37-41]. 

In this paper we examine the 'cosmological radion' which is defined as the inter- 
brane distance of the cosmologically generalised Randall-Sundrum two brane model. 
The radion is in general time dependent and in this paper we go on to investigate this 
time dependent behaviour fully. Many authors impose a stabilisation mechanism on 
the radion such as [42-46], in order to have a static size of extra dimension. While 
this is understandable, most of the mechanisms are inserted in an artificial man- 
ner and it is therefore interesting to examine the nature of the radion without such 
purpose built constraints. Several linearised approaches have been made [45,47-50] 
however, we will define the radion in a non-perturbative way, and then use the as- 
sumptions of homogeneity and isotropy to derive the non-linear equations of motion 
of the cosmological radion, thus providing a more general and elegant proof to that 
found in [51]. 

The radion equations will be used to identify scenarios where stable and unstable 
interbrane distances exist, and then the equations will be linearised to determine the 
nature of these equilibrium positions. We then go on to solve numerically for the 
trajectories of the radion in several cases of interest including when both branes are 
dS or AdS, when the reference brane possesses a sensible cosmology, when both 
branes are either radiation or matter dominated, when phase transitions occur, and 
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when the Z2 symmetry is broken allowing scenarios with two positive tensions branes. 
Hence we significantly extend the results found by [51], in which only fairly trivial 
situations were looked at. A rich variety of behaviour of the radion is found and in 
order to understand this we present a conceptual argument based on the perspective 
of a bulk observer, summarised in section 

2. Deriving the Complete Non-Linear Equations for the Cos- 
mological Radion 

In this section we first describe the generalised Randall-Sundrum two brane scenario 
before defining the cosmological radion in a non-perturbative manner. The complete 
non-linear equation for the cosmological radion or interbrane distance 7?.(t) is then 
derived for the general case. It is assumed that no Z2 symmetry exists across either 
brane and that the bulk possess a non-zero Weyl tensor component and we there- 
fore generalise the derivation presented in [51,52]. Alternative, mainly linearised 
approaches can be found in [45,47-50]. 

2.1 Cosmological Two Brane Scenario 

Since we are interested in cosmological solutions, we take a metric of the form: 

ds^ = -n^{t,y)^t^ + o?{t,y)-iij^x'^x^ + d?/^ (2.1) 

where y is the coordinate of the fifth dimension and we adopt a brane-based approach 
where the reference brane is the hyper-surface defined by y = 0. Here is a 
maximally symmetric 3-dimensional metric with k = —1,0,1 parameterising the 
spatial curvature. The metric is found by solving the five-dimensional Einstein's 
equations which take the form. 

Gab = Rab — -Rqab = ^^^Tab, (2.2) 

where Rab is the five-dimensional Ricci tensor, R = q^^Rab the scalar curvature 
and where we define = 1/M|, M5 being the fundamental (reduced) 5D Planck 
Mass. The stress-energy-momentum tensor can be written as, 

T^B = T^Blbrane + T^slbulk- (2.3) 

Again with cosmology in mind, we assume a homogeneous and isotropic geometry in 
both branes and this makes it possible to write the first term as: 

T^Blbrane = S {y) dieig{- po , Pq , Po , Pq , 0) + (5(y - 7^(^)) diag(-p2, ^2, ^2, A, 0), 

(2.4) 

where po and Pq are the energy density and pressure of the reference brane, with 
similar definitions for the second brane. Here we have defined the position of the 
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second brane to be at y = 7^(t) which is in general time dependent. 7?.(t) is the inter- 
brane distance otherwise known as the cosmological radion, and will be investigated 



throughout this paper. The second term in equation (2.3), which describes the same 



negative bulk cosmological constant A either side of the brane, is of the form, 

T^Blbuik = ^ diag(-A, -A, -A, -A, -A). (2.5) 



Substituting the metric given by equation (|2.1|) into the Einstein equations (|2.2| ) 
allows one to find the metric coefficients a{t, y) and n(t, y) in the bulk [53], which are 
used throughout this paper. The Friedmann equation on the reference brane can also 
be found by using the Einstein equations combined with the Israel junction conditions 
which relate the jump in the extrinsic curvature tensor to the energy-momentum- 
tensor on the brane [54-57]. This gives the well known (non-Z2 symmetric) brane 
world Friedmann equation [15,53] (with = a{t,y = 0)), 



' = ^7o) = 6 + 36^° ~ ^ + ^ + ^^-^^ 



where C and F are constants related to the Weyl tensor and the lack of Z2 symmetry 
respectively. These constants are in fact related to the black hole masses on either 
side of the reference brane by C = (Co + Ci)/2 and F = 3(Ci — Co)/4k^, where we 
have defined the black hole mass to the left of the reference brane {y < 0) as Cq and 
to the right {0 < y <7l) as Ci. Below we will also use the mass of the black hole to 
the right of the second brane {TZ < y) defined as C2. 

In order to avoid the redundancy of several physical constants it is desirable to 
define a mass scale /i and an energy density scale a such that, 



/ = . = %■ (2.7) 



We can then deal solely with the dimensionless energy densities (r/o and 172) and 
pressures (po and P2) of the branes, defined by, 

^0 = — , ^2 = — , (2.8) 
o a 

Po = — , P2 = — • (2.9) 
o a 

These dimensionless quantities will be used throughout the rest of this paper. We 
now go on to derive the non-linear equation for 7^(t). 

2,2 An Elegant Method for Finding the Equation for il. 



Using the dimensionless quantities defined by equations ( P^ and ( p.9[) , the Fried- 
mann equation on the reference brane now looks like, 

iil{t) = Ao^-P^-^ + ^ + §:^, (2.10) 





a, 



2aQ IGfi'^riQal 
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where Cq and Ci are the masses of the Schwarzschild black hole either side of the 
brane. If the second brane follows a trajectory given hj y = TZit), then the induced 
four-dimensional metric on the second brane is given by, 



n\t,n{t)) -n^ dr + a\t,n{t))dx\ (2.11) 

dr2 + a^(r)dx^ (2.12) 



where the dot represents the derivative with respect to the reference brane time t, 
the subscript 2 implies evaluation on the second brane at y = 71, and r has been 
defined as the proper time as seen by an observer on the second brane and is defined 
by, 



dr = yn2(^,7^(^)) -7^2 dt. (2.13) 

The expansion rate of the second brane as seen by an observer on our brane is simply, 

/ ^ 1 dao ( a a' - A 

H2{t) = -— ^ = - + -7^ , 2.14 
02 at \a OL J 2 

where the dash represents the derivative with respect to and as seen by an observer 
on the second brane itself, 

- - (^^'^)^("^-^or- 

This equation for 7^2 is important as it relates the expansion rate of the second 
brane (with respect to its proper time) to TZ and IZ. It does not seem so useful at 
first, as calculating 0,2 could be difficult. However, it is known that the brane world 
Friedmann equation given by (|2.10| ) is derived from a purely local analysis and that 



should we have chosen the second brane to be stationary and at y = 0, we would 
have derived the equivalent Friedmann equation with ?7o replaced by r]2. This means 
that 7^2 (t) must have the following form, 

n_,2( \ 2 2 2 k C1 + C2 {C2-C1Y /^,1'7^ 

where we have defined as before rj2 = P2/0" and Ci and C2 are again the masses of 
the black hole to the left and the right of the second brane. It should be understood 
that equation (|2.17| ) ensures that the second brane evolves according to the junction 



conditions of the extrinsic curvature tensor, just as equation (|2.10|) ensures similar 
behaviour for our brane. At this point we use the following bulk identity which is 
obtained from the Einstein equations and derived in [12], 

— =— + 2.18 
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where as before z = 0, 1 or 2 for the regions y < 0,0 < y < TZ, or y > TZ respectively. 
Combining this with equation ( p.l7| ) gives, for the < y <TZ case, 



7^9 



a 

na 



— + 

2 

2 ^2 



fiV2 + 



f2 



(2.19) 



where the usual notation /2 = (C2 — Ci)/4/i has been used. Substituting ( |2.19| ) into 
( |2.16| ) and rearranging gives the following equation for TZ, 

1/2 



a 

- + 

a 



71^ a 



n 



IJ'V2 + 




V2ai 



(2.20) 



where we have dropped the subscript '2' as from now on all the metric components 
are to be evaluated at y = 71 unless stated otherwise. This therefore generalises 
the result derived by [51]. Equation ( |2.20|) can be algebraically solved to give the 
generalised first order equation for TZ: 



n 



n 



a'a 

a?"n 



± H2 /U?72 + 



/2 



-1 



(2.21) 



^2ai/ J 

One has to be careful when considering in which situations solutions exist and in 
some cases which sign should be chosen in equation (|2.21|) . We need to consider both 
equations (|2.20|) and (|2.16| ), a fact that was neglected by [51] and led to them having 
to change the ± sign in equation (|2.21|) by hand in order to enable their numerical 
computation. Equation ( p.20| ) tells us immediately that —n < IZ < which just 
means that the motion of the second brane relative to the first cannot be greater 
than the local speed of light. With this in mind the same equation now tells us that 
if a' /a and iJ,rj2 + f2/V2<^2 of the same sign, then if, \a'/a\ > \a/na\ there will be 
two solutions, however if |a'/a| < \a/na\ only one solution exists. On the other hand, 
iia'/a and ^j^t]2 + f 2/^2(^2 ^I's of the opposite sign, \a' /a\ > \a/na\ implies no solutions 
and \a'/a\ > \a/na\ implies only one. For the cases where there is one solution 
only, the overall choice of the ± in ( p.21j ) should be the sign of d(yU?72 + f2l'U2<A)- 
Inspecting equation (|2.16|) we see a similar story: if ?-^2 and a/na are the same sign, 
then if |a'/a| is greater/less than \a/na\ it implies one/two solutions. If 7^2 and d/na 
are of opposite signs then if \a' /a\ is greater/less than \d/na\ it implies one/zero 
solutions. The overall sign for the single solutions in this case should be that of 
1-120! I a. Therefore these two equations compliment each other: when (p.20|) implies 
there are two solutions, equation (|2.16|) shows that there is in fact only one (which 
one depends on the sign of 7^2)- This was neglected by [51], and will be used in 
generating the numerical solutions in section 

2.3 Derivation of the TZ Equation 

In order to investigate equilibrium positions of the second brane it is necessary to 
derive the equation for IZ. Taking the time derivative of equation ( |2.20| ) and then 
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adding to it H2 times equation ( p.20|) gives 




where the function g{a,'r]2) is given hj g = fir]2 + /2/^2«'^- Noting first that 



dt 



37 = l^V2 



f2 m _ ^ /2 Ida 
?72a^ a dt ' 



(2.23) 



and the identities: {l/a)da/dt = H2 and rj2 = —3H2{r]2 + P2) which is just the 
energy-momentum conservation equation for the second brane, one then finds that 
by multiplying ( p.22| ) by 1 — 'R? /n? and rearranging, one gets, 

- i7^ (G55 + Goo/^') (1- 




1-2- 



72^ 



H2 -2nr]2 - 3np2 + 



2^4 I 



ma 




(2.24) 



We will be mainly interested in the case where the bulk possesses a cosmological 
constant, and therefore from the Einstein equations ( p.2| ) we have Gqs = and 
Goo/n"^ = —^55 = — This means that the first two terms on the left hand side of 
equation (|2.24|) vanish. In fact a more general argument first proposed by [51] shows 
that these terms will vanish as long as the bulk energy-momentum tensor satisfies 
the three-dimensional symmetries of homogeneity and isotropy. >From (|2.24|) we 
can now see that the generalised equation for TZ is given by, 



n 



4 + 




hTZ 

n 71^ 



3/2 

-2/ir^2 - 3/ip2 + (1-^) . (2.25) 



Equations ( |2.21| ) and ( |2.25| ) govern the evolution of TZ and will be used to examine 
the behaviour of the general cosmological radion. 



3. The Linearised Radion Equation 

In this section the full equations (|2.21|) and (|2.25|) that govern the evolution of the 
cosmological radion are examined analytically. Unfortunately, due to the highly non- 
linear nature of these equations finding an analytic solution is impossible in all but 
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the most trivial of cases. Therefore, here the stable solutions will be identified and 
the nature of their stability discussed along with a few other important features that 
the solutions must possess. These results will then be compared with the extensive 
numerical analysis that occurs in the next section. 

3.1 Equilibrium Solutions 

In order to find the equilibrium solutions, TZ and TZ are both set to zero in equations 
( |2.21|) and (|2.25|) . This restricts the equation of state on the second brane to be of 



a specific form. It will be shown in general that in order for the size of the fifth 
dimension to be constant, an equation of state of the form p2 = inevitably leads 
to a time dependent U2. Setting 7^ = = in equation (|2.21| ) leads to. 



— = f^V2 + (3.1) 

where a and a' are evaluated at position y = TZ and we have defined r/2 to be the 
dimensionless energy density on the second brane at equilibrium. Rearranging gives. 



2/ir^2 = - ± a/— --4^. 3.2 

a V a a* 

This shows that if for example a'/a < for all y (as is the case for a cosmologically 
realistic Z2 symmetric brane where a(t, y) = cosh fJ,y — riQ sinh fxy and ?7o > 1) and the 
non-Z2 symmetry breaking parameter satisfies < /2 < {aa')"^ / ifj, then the possible 
solutions for f/2 are both negative. The case where /2 > {aa'Y/ifi of course has no 
equilibrium solution. If however, /2 < then one of the solutions will be positive, 
therefore one could have two positive tension branes in equilibrium in a semi-infinite 
five dimensional space-time. This corresponds to the setup examined in [58] and will 
be discussed in more detail later. Setting 7?, = 7?. = in equation (|2.25|) leads to a 
similar constraint on p2, 

n' _ _ 3/2P2 

— = -2fi7]2 - 3fip2 + (3.3) 

where p2 is the dimensionless pressure on the second brane at equilibrium. This 
constraint can be rewritten using equation ( |3.1| ) to give. 



f)2 (- + 2/1772) . . 

3 (- - 2fxri2) 

This demonstrates that except for a few trivial situations, the equilibrium require- 
ment forces UJ2 to be time dependent and therefore for the second brane to have an 



'unnatural' equation of state. Note that equations (|3.2|) and (3.4) do reproduce the 



Randall-Sundrum conditions r/o = —i]2 = 1 or equivalently K^po = — /t^P2 = V— 6A, 
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in the constant brane tension (wq = ^2 = ""1), non-Z2 and time independent limit 
as expected. For the specific case where the equation of state of the reference brane 
is given by po = with u;o = —1, one finds from solving the Einstein equations 
that a! I a = n'/n. Equation ( p.4|) then becomes (using equation 




Therefore, at times when the non-Z2 symmetric nature of the second brane dominates 
i.e. when /2/?72a^ ^ fJ^V2 equation (|3.5|) implies that U2 = 1/3. This corresponds 
to late time in a radiation dominated universe. Alternatively, if the second brane is 
approximately Z2 symmetric such that f2/f]2Ci^ <^ yU?72 then equation ( p.5|) predicts 
that UJ2 = —1, which corresponds to the second brane possessing a constant brane 
tension. Another interesting case where equilibrium positions exist is when both 
branes are tuned, Z2-symmetric and possess realistic cosmologies. Here, equilibrium 
positions exist provided both branes have the same equation of state. This situation 
is discussed fully in section ^]4|. The nature of the stability of the above equilibrium 
solutions will be examined in the next section. 



3.2 Radion Fluctuations 

In order to examine the nature of the equilibrium solutions found above it is necessary 
to linearise the radion equations ( p.21| ) and ( p.25| ) around the equilibrium points. 
Setting TZ = STZ, IZ = 6TZ, rj2 = fj2 + ^^2 and P2 = P2 + Sp2 results in, 

- + = - ^^2 - —^oa , 3.6 

a J an n \ Vi^^ J ^20 

n J 77,^ n ri2a^ \ P2 ^2 ci J 



In order to evaluate the first term in each of the equations ( |3.6|) and (|3.7|) it is now 
assumed that the reference brane respects Z2 symmetry and possesses no Weyl tensor 
component. Therefore /o = and C = or in terms of the bulk Schwarzschild masses 
on either side of the reference brane: Cq = Ci = 0. Since C2 7^ the second brane 
will not in general be Z2 symmetric. This means that the metric components, which 



are obtained by solving Einstein's equations (|2^ ), now look like [51], 



a{t,y) = ao(cosh/i?/ — ?7o sinh/iy) , (3.8) 
n(t,y) = cosh fiy — fjo sinh fiy , (3.9) 
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where tjq = rjo + fjo/ Hq. It can now be explicitly shown using equations (|3.1|) , (|3.3| ) 

/2 



(Ml) and (U) that, 
(5 



6 



/i - -2/ir72 



- 3/ip2 



57^ 



ra„5TZ . 



57^ = m^57^, 



and defining the function g(t,7V) as, 

^(t,7^) 



/2 



^20"* 

equations ( |3.6| ) and (|3?7|) can be rewritten in the following form: 

n2— + {ml + 4g{g~fif]2))Sn = (2/i - ^)5r/2 , 
n r/2 

— - — — + {ml + 12^g{g - fir]2))Sn = 



(3.10) 
(3.11) 

(3.12) 

(3.13) 



n 



V2 



,P2. 



g 



-2(/i + 3^{g - fif]2))6v2 - 3(2/i - ^)5p2 ,(3.14) 

V2 m 

where the fact that Sa/a = {a' /a)51Z and that {a/an)6'JZ = T-C2S'JZ to first order, 
has been used. Thus, it can be seen that in general the radion fluctuations are 
entangled with the matter fluctuations and that solving the linearised equations 
becomes non-trivial. It is possible however, when the fluctuations are of the adiabatic 
type, 6p2 = c\5r]2 to take the appropriate linear combination of equations ( p.l3| ) 
and ( p.l4|) to obtain an equation that depends upon the radion fluctuations only. 
Examining initially the case where ujq = —l,p2 = UJ2V2 and therefore = 002, 
and where the second brane is dominated by its lack of Z2 symmetry such that 
12/ '^20''^ ^ At'72, ones finds that equations (|3.13D and (|3.14| ) combine to give. 



571 rr — 3002^2^71 . 



f2 

+ 3u;2(l - 30^2)^^ 



6n = 0, 



(3.15) 



where only the dominant terms have been kept, r is defined as before as the time 
experienced by an observer on the second brane. The matter fiuctuations are now 
related to the radion fiuctuations by the formula: 



Sri2 



f2 



7i26nr + 



(3.16) 



Equation ( |3.15| ) can now be solved for a radiation dominated second brane where 
UJ2 = 1/3. Using the approximation 7-^2 — —f 2/^2(1'^ and remembering that in this 
situation /2 < 0, one obtains the expression, 

,4 



6n 



^2i/2 



(3.17) 



10 



where the subscript i means the value of the quantity at Xj. Therefore it is clear that 
when there is a constant brane tension on the reference brane and a highly non-Z2 
symmetric second brane such as is the case suggested by Randall and Lyken [58], 
an equilibrium solution exists only when the second brane is radiation dominated, 
however this solution for the interbrane distance or cosmological radion is inherently 
unstable. 

Another case of interest is when both branes are Z2 symmetric and as discussed 
in the previous section if cuq = —1 then an equilibrium solution exists only when 
U2 = —I also. Adding three times equation ( p.l3|) to equation ( ^.14 ) and neglecting 
all /2 terms results in, 

^2 

6n,rr + 3H257^,^ + Afi'^{l-r]'^)6n = -y^^2, (3.18) 

where ST2 is the linear perturbation of the trace of the matter energy-momentum 
tensor [51] and is given by, 

—ST2 = ~ij6r]2 + 3/i5p2- (3.19) 

It can then be seen that equation (|3.18|) is just the equation of motion of a scalar 
field coupled to the perturbed trace of the matter energy-momentum tensor. This 
therefore confirms previous results found in [47-49]. The coefficient of the STZ term 
in equation ( p.l8|) demonstrates that when both branes are (iS'4 implying that both 



?7o > 1 and 773 > 1, small perturbations around the equilibrium radion position are 
unstable. Alternatively, if both branes are AdS^ and therefore r^Q < 1 and t]"^ < I 
then the perturbations are stable. It is again possible in this case to eliminate the 
matter fiuctuation terms by taking a suitable linear combination of equations ( p.l3| ) 
and ( p.l4|) . For a general U2 and keeping uq = —I this leads to, 



SUrr + {2 + 3u2)n2 Sn^r + 3/1^(1 + c^s) (1 " (2 + 3tU2)r/2) 57^ = 0. (3.20) 

The matter fiuctuations are now related to the radion fiuctuations by, 

I^Sr]2 = 7^2 + /i2(l -r/^) 57^. (3.21) 

Assuming that both branes are dS^, that uj2 = —I and that k = Ci = C2 = ensures 
that TC2 is time independent and given by l-i2 = l^^ijll ~ It is therefore possible 
to solve equation ( p.20[ ) for TZij) giving. 



57^ = 57^,e^^(^-^") - -^57^2,(1 -e«^(^-^')). (3.22) 

/T2 



where as before SlZi and 5rj2i are the initial radion fluctuations and energy density 
fluctuations at time r^. This describes the behaviour near the equilibrium position 
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and demonstrates explicitly the unstable nature of the radion for two dS4 branes and 
will be compared with several numerical solutions discussed in the next section [51]. 
If instead, both branes are assumed to be AdS4^ then it is not possible for H2 to be 
time independent and either k, Ci or C2 must be non-zero. This leads to a recollapsing 
universe on the second brane. Again taking 002 = —I allows equation (p.20|) to be 
written as simply, 

Sn^rr = ^2 SUr- (3.23) 

Examining equations ( p.21j ) and ( |3.23| ) it can be seen that the radion will only exhibit 
a pseudo-stable behaviour in that TZ will initially accelerate toward the equilibrium 
position but will only decelerate once H2 has changed from positive to negative. 
In the section |5| the analytic results discussed here will be compared with complete 
numerical solutions of the radion equation for several cases of interest. Before this we 
go on in the next section to discuss the cosmological radion in terms of the static bulk 
perspective, in order to obtain an intuitive understanding of some of the numerical 
results that are presented in section ||. 



4. The Two Brane Scenario: a Bulk Perspective 

In the previous section we linearised the radion equations in order to determine the 
stability of several equilibrium positions, and therefore to obtain some idea of the 
behaviour of the radion in several cases of interest. Here we go on to briefly examine 
the two brane model from the perspective of a bulk observer. This will allow us to 
show that the evolution of the radion can be thought of as a competition between 
the expansion rates of the two branes, and this will help determine in which cases 
the branes collide or move apart. These predictions will be shown to be in agreement 
with the numerical results presented in the next section. 

We are considering purely five-dimensional spacetimes that respect the sym- 
metries of three-dimensional homogeneity and isotropy, and that only possess a five 
dimensional cosmological constant in the bulk. It was first shown by [15] that a modi- 
fied version of Birkhoff 's theorem implies that such spacetimes are in fact static. The 
bulk metric takes the general form: 

ds^ = -(f){R)dT^ + R^-fijdx'dx^ + 0(i?)-Mi?^ (4.1) 

where we have defined the function 4>{R) as, 

0(i?) = (^f^'R' + k - (4.2) 

and as before k takes the values 0, —1, or 1 for fiat, closed or open geometries, and 
7jj is the corresponding metric on the unit plane, hyperboloid or sphere. T is the 
time experienced by a bulk observer, and R the coordinate of the fifth dimension. 
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while as before the subscript i is either 0, 1 or 2, corresponding to the bulk regions 
to the left, between or to the right of the branes. 

The simple form of these bulk solutions is counter-balanced by the fact that in 
these coordinates both the 3-branes cannot (except in trivial cases) be stationary, 
and in general move through the bulk spacetime. In fact in the non-Z2 symmetric 
case, the brane moves through two different bulk spacetimes, with for example, two 
different bulk masses either side of the brane. If the reference brane has a trajectory 
given by -R = Rq{T), then the induced metric on it is given by, 

ds^ = -dt^ + Rl{t)-fijdx'dx^, (4.3) 

where again t is the time experienced by an observer on the reference brane and is 
related to the bulk time T by, 



df 



0(i?o) - 



1 fdRo^^ 



0(i?o) V dT 



dT\ (4.4) 



Hence from equation (|4.3| ) one can see that the apparent expansion of the brane 
as seen by such an observer is in fact caused by the brane's motion through the 
bulk [15]. One can therefore deduce that the brane's position -Ro('^) in these bulk 
coordinates is proportional to the scale factor ao{t) of the brane used throughout the 
rest of this paper. Thus the trajectory of the brane, at least in terms of t, can be 
found from its Friedmann equation ( p.lOj ). In order to compare the trajectories of 
two such branes we need to solve for Rq in terms of the bulk time T. Equation (|4.4| ) 
can be rearranged to give, 

d^ ^ 0(i?o) . . 

dT {<p{Ro) + H^,Riy/^' ^^-"^ 

where Hq is again the Hubble parameter on the reference brane given by HqRq = 
dRo/dt. This can be combined with equation ( |4.5| ) to show that the brane's motion 
in terms of T satisfies, 

dRo ^ dRo d^ ^ HpRo^iRo) 

dT dt dT {<P{Rp) + HlRiyl'^- ^ ■ 

By replacing Rq and Hq with R2 and if 2, one can obtain the equivalent equation for 
the second brane, 

dR2 ^ n2R2<P{R2) , . 

dT {<P{R2)+'HlRiyi^' ^ 

It must be remembered that here Rq{T) and R2(T) correspond to the positions of 
the branes in the static spacetime that lies between them. The interbrane distance 
in terms of these bulk coordinates is now just given hj R = Rq — R2, where Ro(T) 
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and R2{T) can be found by solving equations (|4.6| ) and (|4.7|) . For example, if we 
assume the reference brane is Z2 symmetric, equation ( |4.6| ) becomes, 



Ho4>{Ro 



dT /i|?7o| 

which assuming that k = Cq = Ci = 0, reduces to. 



(4.8) 



dRo 
dT 



(4.9) 



If the reference brane possesses a constant dimensionless energy density rjo, then this 
can be solved to give the position of the reference brane to be. 



RoiT) 



1 

Roi 




-1 



(4.10) 



where Roi is the position of the brane at time Tj. If we make similar assumptions for 
the second brane, in that it is also Z2 symmetric (C2 = 0) and that 772 is constant, 
then R2(T) is found to be of the same form. 



R2{T) 



1 

R2i 




-1 



(4.11) 



where again R2i = R2{Ti). Now taking the convention i?oi > R2i, which here cor- 
responds to the reference brane being of positive tension and the second brane 
being of negative tension, we can ask under what conditions the branes collide. 
Ro(T) = R2(T) occurs at a time Tc given by, 



^/R2i — V-Rpi 

l^{T-C2/\V2\ - Hq/tIq 



(4.12) 



however the collision must occur for positive Rq and R2, and Ro(T) > only when 
T < Too, where Too is given by. 



T 



^0 



^HoRoi 

Therefore a collision will occur only if Tc < Too, which implies the condition, 



(4.13) 



R2i 
Roi 



> 



Ho\r]2\ 



1 - 1/77^ 



1/2 



(4.14) 



So using the bulk perspective we have found the requirement for brane collision to 
occur in terms of the branes' energy densities and initial positions, in the constant 



brane tension case. This condition ( |4.14| ), is in fact a specific case of the brane based 
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equilibrium conditions and (|3.3|) . This will be discussed further in section |5]T 



where we look at the constant brane tension case from the brane perspective, and 
where the equivalence of these conditions will be demonstrated explicitly. 

For more realistic cases such as where both branes possess a cosmologically evolv- 
ing energy density, the bulk perspective is less useful. It is, in general not possible 
to analytically solve equations ( |4.6|) and ( ^.7| ) in order to obtain the trajectories of 
the branes and hence it is difficult to determine whether the branes will collide. In 
addition one cannot easily determine the conditions for equilibrium of the interbrane 
distance; however it must be noted that this equilibrium is in some sense a creation of 
the brane coordinates themselves. One inevitably has to turn to numerical methods, 
and due to this in the next section we numerically integrate the radion equation, 
formulated in terms of the brane coordinates. Note that in the brane based perspec- 
tive, all the information about the branes' relative positions contained in equations 
( [4.6| ) and ( ^.7| ) is described by one equation ( p.21| ), the non-linear equation for the 
radion. Hence we only need to numerically solve one non-linear equation. 

The bulk perspective however, can still provide intuitive insight into the results 
obtained in the following chapter. Thinking of the evolution of the interbrane dis- 
tance as a competition between the branes' expansion rates helps explain the rich 
variety of behaviour that is found. 



5. Numerical Analysis of the Non-Linear Radion 

Previously, the equations for the cosmological radion have been derived and anal- 
ysed. The equilibrium points have been identified and the stability of these solutions 
examined. Unfortunately, due to the highly non-linear nature of the radion equation 
there is only so much one can learn from analytical methods. Hence, in this sec- 
tion, the radion equations are solved numerically in several different cases in order 
to determine the non-perturbative aspects of the cosmological radion. The subtle 
numerical problems involved are discussed, and the results are interpreted using the 
bulk perspective described in the previous section. We first review and then extend 
some of the results on dS and AdS branes presented in [51]. 

5.1 Constant Tensions on Both Branes: The dS Case 

If the assumptions are made that the reference brane is dS^ or equivalently that 
r^o > 1, and that Co = Ci = implying that the brane possesses Z2 symmetry and 
has no Weyl tensor component, then the brane based metric components can be 
written in the simple form, 

a{t,y) = ao{t) {cosh fiy - r]osmh^y) = ao{t)^r]^ - Isinh fi{yh - \y\) , (5.1) 
n{y) = {cosh fiy -r]osmh fiy) = J r]"^ - 1 sinh fi{yh - \y\) , (5.2) 
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where yu denotes the position of the coordinate singularity defined by a(t, = and 
is given by, 

avh = tanh~^ — . (5.3) 

^0 



Replacing the metric components (|5.1|) and ( |5.2| ) into equations (^]T|) and ( |3.3| ) to 
determine the equilibrium radion position denoted by j/e, one finds as stated before 
that the equation of state on the second brane must be of the form p2 = —^2 and 
that, 

UDe = tanh~^ — + tanh"^ — . (5.4) 

In order to numerically solve the radion equation here and in the following sections, 
it is useful to convert to the following dimensionless variables, 

Z = /i7^, Zh^e = fJ'yh,e s = fit, ho = — , /12 = — . (5.5) 



Replacing the expressions for the metric coefficients ( ^.1|) and (|5.2|) into the radion 
equation ( p.21| ) and changing to the above dimensionless variables gives. 



, _ u . cosh(2:;, -z) + h2V2 smh\zh - z) 

where it should be remembered that here = r^Q — 1, and setting k = gives 
^2 = ^2 ~ 1- Since the energy densities of both branes are by assumption constant 
in time, the evolution of the radion is completely embodied by equation (|5.6| ) and 
all that is required to solve it is the initial radion position denoted by Zi. Numerical 
solutions were generated for two different situations, each with varying initial radion 
positions and the results are shown in figures |l] and |2|. Figure || corresponds to the 
case where ?7o = 1.5 and 772 = —3 and shows the trajectories of six different initial 
radion positions distributed around the equilibrium position Zf,. Equations (|5.3| ) 
and ( p.4|) show that for this choice of brane tensions the equilibrium and singularity 
positions are Ze — 0.458 and Zh — 0.805 which correspond to the first and second 
dotted lines of figure |l]. As expected from the above analysis, the equilibrium position 
is found to be unstable: initial positions satisfying Zi < z^ lead to the second brane 
colliding with the reference brane; as opposed to initial positions satisfying Zi > z^ 
when the second brane asymptotically 'freezes out' at Zh as seen by an observer on 
the reference brane. Consideration of equations ( ^.61) and (|2.13|) demonstrates that 
crossing z^ only takes a finite amount of proper time on the second brane. 

The situation where rjo = 1.5 and 772 = —1-3 was investigated and the results are 
shown in figure Referring to the expression for the equilibrium position given by 
equation (|5.4| ), it can be seen that if both branes are dS with for example 770 positive 
and ?72 negative, then |?7o| > \ri2\ implies that Ze < 0. Due to the Z2 symmetry this 
means that there is no physical equilibrium position for the radion in this setup. This 
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Figure 1: Both branes are dS and tjq = 1.5, rj2 = —3. 
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Figure 2: Both branes are dS and r/o = 1.5, r]2 = —1.3. 

is confirmed by figure § where every choice of initial position for the radion leads to 
the second brane freezing out at Zh, where Zh — 0.805 as in the previous example. 
This behaviour can be better understood from a bulk observer's perspective. In the 
bulk coordinate system, both branes are moving through a static background. Their 
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trajectories from this point of view are solely dependent upon their expansion rates 
or Hubble constants; increasing the magnitude of 772 increases 0,2 and hence improves 
the chances of the second brane 'catching' the reference brane. Alternatively, if as 
in figure 0, the magnitude of 770 is greater than that of 772, it will be impossible 
for the branes to collide. In section || we used the bulk perspective to solve for the 
trajectories of two constant tension dS branes, and derived the condition for collision 
( [4.14|) : R2i/Roi > HQ\r)2\/H2f]o- This can be rewritten using the fact that Rq and R2 
are proportional to ao{t) and a{t,z) respectively, which implies, 

aa{t,z) Ho\r]2\ 
ao{t) n2Vo' ^ ' 

where we have scaled the three spatial dimensions so that Rq = oq and there- 
fore R2 = aa{t, z) where a is some positive constant discussed below. Again using 
the explicit expression for the metric element a{t,z) given by equation ( |5.1|) , this 
becomes, 

|?72| 1 

smh[zh-z) > — — -. 5.8 

ario {ril - 

Now assuming that a = |r72|/^o, one finds that the branes will collide if the initial z 
satisfies, 

tanh(z/i — z) > - — -, (5.9) 

ml 



from which one can recover the equilibrium condition equation (|5.4| ), derived from 
the brane perspective. One can check that a does indeed take the required form, as 
it must in order for the two perspectives to agree, by checking the explicit transfor- 
mation between bulk and brane coordinates as detailed in [59]. 

The bulk perspective is very useful for understanding the behaviour of the non- 
perturbative cosmological radion and will be discussed further. It must be stressed 
that a time independent bulk solution will not be possible for scenarios that contain 
more general matter in the bulk, such as for example scalar fields or dilaton fields. 
This is another reason for using brane based coordinates to examine the radion. 

5.2 Constant Tensions on Both Branes: The AdS Case 

In this section the radion dynamics are investigated for the case of two AdS^ branes, 
that is two branes such that |?7o| < 1 and |?72| < 1- Again it is assumed that both 
branes are Z2-symmetric and that they possess no Weyl tensor component such that 
Cq = Ci = C2 = 0. This implies that one must take k = —1 to get a consistent 
solution to the Einstein equations. The metric coefficients take the form. 



a{t,y) = ao{t) {cosh fiy - rjosinh fiy) = ao(t) W 1 - r^g cosh (|z| - 2;^) , (5.10) 



n{y) = {cosh fiy - r]osinh fiy) = \ 1 - r]^ cosh fi{\z\ - Zm) , (5.11) 
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where the dimensionless coordinate z = ny has been used as before, and z^ which 
corresponds to the minimum of a(s, z) and n(s, z) is given by, 

Zm = tanh""^ r^o- (5-12) 



From equation ( |2.20|) it is seen that since a(s, z) possesses a minimum, it is possible 



to have two positive tension branes in a compactified five-dimensional spacetime [60]. 
The radion equation ( |2.21|) now looks like, 

z = Kq cosh.{z — Zm) X (5.13) 



hosmh.{z — Zm) ± 1]2 cosh(z — Zm)y + Kq{1 — cosh {z — Zm)) 

hi + rilKl cosli^ {z - Zm) '(5-14) 



where Kq = ^Jl — r]l and K2 = ^/T^-rj^. Using equation ( |2.10| ) the dimensionless 
Hubble parameter on the reference brane is now given by, 



ho = -K^ + (5.15) 

Uq 

which can be solved to give ao(t) = sm{Kos) / fxKo and therefore ho = Ko/ taii{Kos) 
which is all that is needed in order to numerically solve equation (|5.13|) . It must 



be noted that, when generating a solution to equation ( |5.13| ) one must be careful 



to ensure that the sign of the square root is changed whenever the argument of the 
square root vanishes i.e. when /i2 = 0. This has to be the case and can be seen on 
examination of equation ( p. 211) . The equilibrium position can now be shown to be. 



z 



■e 



= tanh ^ r/o + tanh ^ r]2- (5.16) 



Figures ^ and ^ show the radion trajectories that were numerically generated for 
two distinct situations. The first case that was examined, corresponding to figure 
1^, was where 770 = 0.6 and rj2 = 0.6 which implies choosing the positive root in 
equation (|5.13| ). The initial radion positions are distributed about the equilibrium 



point Ze = 1.39 which is represented by the dotted line. As expected from the analysis 
in section |3]^, it is seen that the the radion displays a pseudo-stable behaviour - 
initially the interbrane distance accelerates toward z^, but it only begins to decelerate 
once /12 = which does not necessarily correspond to z crossing the equilibrium 
point. It is this fact combined with the warped nature of ?t.(s, z) that causes different 
trajectories to reach /12 = at different points and hence leads to the asymmetrical 
(in the z direction) appearance of figure ||. Figure | shows several radion trajectories 
when the second brane possesses a negative tension. Specifically, the brane tensions 
have been chosen to be rjo = 0.9 and 772 = —0.9 which implies that Ze = 0. As 
expected, this leads to the second brane colliding with the first for any given initial 
radion position. 
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Figure 3: Both branes are AdS and 770 = 0.6, 772 = 0.6. 




Figure 4: Both branes are AdS and ijq = 0.9, r)2 = —0.9. 



5.3 Radiation and Matter on the Reference Brane 



We now turn to a more interesting situation where there exists a realistic cosmology 
on the reference brane. We assume that the dimensionless energy density 770 can be 
decomposed into two parts [61]: a constant brane tension and a 'physical' energy 
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density denoted by r^op such that, 



?7o = 1 + r]op, (5.17) 

and similarly the dimensionless pressure is decomposed into, 

Po = -1 + Pop. (5.18) 

This is done as before to ensure that the reference brane undergoes standard cosmol- 
ogy at late times or more specifically that the Friedmann equation has the form, 

^0 = Ao'p + 2/i2r/op, (5.19) 

so that Hq oc rjop at late times, where we have assumed that both branes are Z2 
symmetric, possess no Weyl tensor component and that k = 0. If the reference brane 
now has an equation of state pop = uJoVop, equations ( p.8|) and (|3.9|) show that the 
metric coefficients in terms of s and z will now be given by, 

a{s, z) = ao(cosh z — {1 + rjop) sinh z) , (5.20) 
n(s, z) = coshz — (1 — (2 + 3co'o))?7op sinh z , (5.21) 

where now both oq and ?7op are functions of the dimensionless time parameter s. 
The Friedmann equation ( ^.19| ) can now be solved for rjop using the brane energy 



conservation equation, 

r]Q = -3HQ{r]o + Po), (5.22) 

which leads to, 

^op(s) = 12 2, ' (^-23) 

where we have defined go = 3(l + ci;o)- Using this expression for riQp{s), the position of 
the coordinate singularity defined by a(s, Zh) = is now given by the simple formula. 



Zh 



= tanh-M-^ — ) = ln(goS + l). (5.24) 



We now assume that the reference brane is in a radiation dominated phase in that 
Wo = 1/3, go = 4 and therefore r]op oc Oq"^ and that the second brane has a constant 
brane tension. It is now possible, using equations (|5.17|) , ( |5.19|) , ( ^.20| ), (|5.21|) and 



( |5.23|) to numerically solve the radion equation (p.21|) . This was done using several 



different initial conditions and the results are shown in figure ^ along with the time 
dependent position of Zh denoted by the dotted line and given by equation ( |5.24| ). 



The six initial radion positions were distributed evenly between z = and z = Zh- 
In figure || the radion positions are graphed against time for a second brane tension 
?72 = —1.2, and as can be seen, all the radion trajectories lead to the two branes 
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colliding. Trajectories that initially begin closer to take longer to return to the 
reference brane as would be expected. 

Similar results were also obtained for various values of the second brane tension. 
It was found that the radion trajectories exhibit the same general behaviour, but the 
second brane would return and collide with the reference brane sooner/later if the 
magnitude of 772 was greater/lesser. An interesting case is shown in figure |, where 
the reference brane is again radiation dominated, however the second brane is now 
a tuned or critical brane with r/o = —1. The initial radion positions were spread 
between z = and z = 0.95zh at s = 0.05 and the coordinate singularity moves 
away from the reference brane as before with Zh = ln(4s + 1). In this situation the 
second brane does not return and collide with the reference brane, instead the inter- 
brane distance increases with logarithmic behaviour similar to that of the coordinate 
singularity at Zh. If instead, it is assumed that the reference brane is in a matter 
dominated phase so that uq = 0, qq = 3 and therefore r^op oc Oq^, the dynamics of 
the interbrane distance will again be altered. The coordinate singularity is now at 
Zh = ln(3s + 1), and it was found that if the reference brane is matter dominated, 
the second brane takes longer to return. 

This general radion behaviour when the reference brane has a more realistic 
cosmology can be understood in terms of the evolution of the brane energy densities 
and by viewing the situation from a bulk observers perspective. Initially rjop is very 
large causing oq to increase rapidly (oq oc t^/"^) and therefore the reference brane will 
be moving rapidly through the static bulk. The second brane on the other hand, 
only has a low constant energy density 772, and will begin by moving slowly, hence 
the interbrane distance will initially increase. At late times, r/op tends to zero and 
the reference brane is much slower compared to the second brane which now has 
a2 oc exp H2T and therefore the interbrane distance will decrease and eventually the 
branes will collide. This means that the time taken to collide is decreased by having 
a larger 772, in agreement with the above results. Alternatively, a matter dominated 
reference brane will maintain a higher 'speed' at late times and therefore will increase 
the time taken to collide. In figure ^, where the tuned second brane never returns, 
we have that ?-^2 = 0, and therefore from the bulk observers point of view the 
second brane is stationary. Since our brane is expanding (albeit with a continually 
decreasing rate), it moves away from the second brane and hence the interbrane 
distance increases monotonically. All these features of the radion's behaviour will be 
discussed further in the next section when we go on to investigate the time taken for 
the branes to collide. 

5.3.1 The Time Taken for the Branes to Collide 

In order to obtain an intuitive understanding of the nature of the radion it is useful 
to investigate explicitly Scoi, the dimensionless time taken for the branes to collide. 
Scoi depends on many variables including the initial conditions, the equations of 



22 



5 10 15 20 

s 

Figure 5: Radiation on the reference brane, 772 = —1.2 on the second brane. 
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state on either brane, the curvature and non-Z2 symmetry of the branes and the 
proximity of the second brane to Zh- Here we have restricted our analysis to the 
situations presented in the previous section, where the reference brane possesses a 
realistic cosmology and the second brane has a constant negative brane tension. 
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The dependence of the collision time on the initial radion position Zj, when the 
reference brane is either radiation or matter dominated was numerically determined 
and the results are shown in figure |^. Firstly one can see that for all values of Zi/zh, 
the branes collide much sooner for a radiation as opposed to a matter dominated 
reference brane. As discussed in the previous section, this is due to the matter 
dominated brane moving through the bulk with a greater 'speed' than that of the 
radiation dominated brane. Secondly, as the initial radion position increases, the 
collision time also increases, becoming singular as Zi — > Zh. The reason for this is 
that a(s, z) — at Zh, and therefore the initial scale factor of a brane close to Zh 
tends to zero also. Since the second brane has to expand from a{s,Zi) to ao{scoi) in 
order for collision to occur, a vanishing a(s, Zi) implies Scoi oo. 

The effect on the collision time of varying the tension of the second brane was 
investigated, and the results are shown in figure |. A similar effect to figure |^ is 
seen in that radiation domination implies faster collision. In addition, it can be seen 
that increasing the magnitude of 772 leads to the second brane moving faster, which 
in turn decreases Scoi- As can be seen from figure when 772 — — 1 then Scoi — * 00, 
since the second brane is effectively stationary from the bulk point of view. 

5.3.2 The Ekpyrotic and Cyclic Models 

The above discussion of the collision times of the two branes highlights an important 
feature of the well known ekpyrotic and cyclic models. In the Ekpyrotic universe [62- 
67] it was proposed that the supposed beginning of our four-dimensional universe is 
in fact the result of a collision between a boundary and a bulk brane. Some of the 
energy of the collision would be converted into a hot big bang, and it was shown 
how one therefore no longer needs to incorporate infiation to explain the current 
observational data. The cyclic model [68,69] went one step further: it was suggested 
that the big bang was due to the collision of two boundary branes, and that these 
branes would collide periodically every several trillion years. This would also explain 
the apparent accelerated expansion of the universe we observe today. 

Both these models, however, needed to employ a potential between the branes in 
order for the brane collisions to occur at the correct times. Although these potentials 
were no more complex than those employed in the most recent infiationary theories, 
they were to say the least ad hoc. We have investigated throughout this paper 
the natural behaviour of cosmologically realistic branes moving in a five-dimensional 
AdS bulk, and as can be seen from the above results, we find that brane collisions 
(at least in this model) are commonplace. We also however, find that the collision 
time will not be significantly greater than the inverse five-dimensional Plank mass^ 
(~ 30Te\^ or above), and therefore to delay brane collisions to the extent proposed in 
the Ekpyrotic and Cyclic models, a certain amount of fine tuning must be involved. 

fact which can be inferred by dimensional analysis. 
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Figure 7: Time taken for branes to collide for varying initial z 



Taking this point further, in general one might expect a five-dimensional model of 
this kind to contain a large number of branes, each containing different equations of 
state. These v^^ould necessarily be colliding over timescales of around Mg"^, so in the 
case of a Cyclic universe the main question to answer would be why our brane has 
not hit anything else in the last thirteen billion years. 

5.4 Radiation and Matter on Either Brane 

We now turn to the interesting case where both branes are Z2 symmetric and have 
time dependent energy densities ?7o(t) and r^2(^)- This allows us to investigate the 
behaviour of the radion when both branes are either radiation or matter dominated. 
There are, however, some subtleties as defining a time dependent energy density on 
the second brane leads in most cases to — > 0) = — 00, which could lead to some 
phenomenological problems. 

Assuming that both brane energy densities can be decomposed into a constant 
brane tension component and a time dependent 'physical' component we have, 

770 = 1 + r^op, Po = -1 + Pop, (5.25) 
for the reference brane, and, 

?72 = -1 - ^72^, P2 = 1 - P2p, (5.26) 

for the second brane where 772^ and p2p are positive quantities. We have assumed 
that Co — Ci — C2 — k — 0, and we want to obtain the Randall-Sundrum two brane 
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Magnitude of second brane's tension 
Figure 8: Time taken for branes to collide for varying r]2 



model at late times, therefore the assumed form of 772 and p2 given by equation ( ^.261 ) 
is the only possible choice. 

The junction conditions across the second brane give the energy density conser- 
vation equation: 

f]2 = -3H2 {r]2 + P2), (5.27) 
which in terms of rj2p takes the same form, 

r]2p = -3H2 {rj2p + P2p)- (5.28) 

If the second brane has an equation of state of the form p2p = uj2'n2p then we can use 
equation ( ^.281) to write, 

where we have defined g2 = 3(1 + ^^2). A similar argument for the reference brane 
gives also, 

which was implicitly used in the previous two sections. 70 and 72 are simply constants 
relating the energy densities and scale factors on either brane. 

One can now ask whether an equilibrium solution for the interbrane distance 
exists in this situation. Using the conditions for equilibrium given by equations (|3T 
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and (|3.3| ) we can write the equation for the possible equilibrium point z^. in the same 
form as equation (^.41), 



tanh ^ 



1 



+ tanh ^ 



1 



;i + r/op) ' (-l-r?2p) 
Replacing ?7op and ?72p using equations (|5.29|) and ( |5.30|) , and rearranging gives 



(5.31) 



1 



1 



o(?o-g2)lnao - -g2ln — 



(22 



ao 



-1 1^2' 



(5.32) 



This shows that in general the required expression for is time dependent (since 
both ao and 02 depend on t) and that therefore no equilibrium position exists. If, 
however, we examine the late time behaviour of this equation then we can use the 
fact that 02/00 = cosh — ?7o sinh ~ exp (— Ze) to first order^. Equation ( ^.32| ) can 
now be arranged to give, 



?2 



12 - go) In ao + ^^\^ 



(5.33) 



This demonstrates that there will be an equilibrium solution at late times if and 
only if go = 92 i-e. if both branes have the same equation of state. Therefore if both 
branes are radiation dominated, go = 52 = 4 and z^. = 2 In (70/72), however, if they 
are matter dominated then go = Q'2 = 3 and = 3 In (70/72)- 

We now investigate the stability of such solutions by examining equation ( |3.13| ), 
which in dimensionless variables and with /2 = looks like. 



^2 Sz^s2 + (1 - 112) ^-2 



(5.34) 



where we have defined S2 = yur, the dimensionless time experienced by an observer 
on the second brane. Using the fact that we have set 772 = —1 — Vip = — 1 — (72/(2)'^ 
gives, 

'72\92 a' 



8^2 



<l2 



5z 



q2 mp (-1 - '72p) 



(5.35) 



where we have used the fact that at equilibrium a' /a = 772- Combining equations 
( |5.34|) and (|5.35| ) and expressing everything in terms of fi2p leads to. 



6z 



,S2 



V2p 



V2p + 2 



1/2 



?2 - 2 + (g2 - l)r]2p) Sz. 



(5.36) 



This shows that the fluctuations around Ze will be stable as long as g2 > 2. The 
explicit expression for ri2p in terms of the time experienced by an observer on the 



^Technically we have had to assume both that t^l and that Ze is not significantly close to the 
coordinate singularity at Zh = tanh^^(l/77o). 
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second brane S2 can be found from the Friedmann equation ( |2.17|) combined with 
equation ( ^.291 ) and takes the same form as equation (|5.23| ), 

+ <l2S2 

Inserting this into the fluctuation equation allows us to solve for Sz{s2) giving, 

^^^^^^ = ^ ; ^M(g2-1)M ' ^^-2^) 

[q2S2iq2S2 + 2)f' 

where A is an integration constant. Therefore at late times when S2 ^ 1, the 
fluctuations will behave like Sz ~ (52-32)*^^^'^^"^'' giving Szr ~ {q2S2)^^^'^ for radiation 
and Szm ~ (^'252)"^^^ for matter demonstrating that the interbrane distance will tend 
to Ze at late times. Incidentally, the early time behaviour is Szr ~ {q2S2)~^^^ and 
Szm ~ iq2S2) '^^^, and although we have already shown that the equilibrium position 
Ze is not constant at this time, provided it moves only slowly relative to the evolution 
of the fluctuations then these early time results will still be of significance. 

Using equations (|5.20D , (p. 21]) , ( p.26| ), and (|5.29|) the radion equation (|2.21|) was 
solved numerically for four different cases, with both branes assumed to be in a state 
of either radiation or matter dominance and the results are shown in figures |10|, [11 



and |12|. In each case the six initial radion positions were distributed evenly between 
and Zh, and the integration was started at s = 0.1. Figure || shows the evolution of 
the radion z against the dimensionless time variable s when both branes are radiation 
dominated in that go = 92 = 4. The proportionality constants have been chosen to 
be 7o = 1.04 and 72 = 1 and therefore the late time equilibrium position Ze given by 
equation ( |5.33|) is Ze — 0.0784 which is represented by the dotted line. At late time 
the stable behaviour predicted above is evident as all trajectories tend slowly towards 
Ze- The expected early time behaviour is also confirmed as initially the trajectories 
rapidly converge, however some of them possess stationary points which implies a 
time dependent solution for Zg. In fact, the minimum of a trajectory corresponds to 
where the time dependent solution for Ze actually crosses the trajectory. Figure p!0| 
shows the radion's behaviour when both branes are matter dominated. Here the 
constants are 70 = 1.02 and 72 = 1 and therefore Ze — 0.0594. A stable behaviour 
is again observed as expected from the above analysis, and the similarities with 
figure ^ are substantial. Note how the rate at which trajectories converge at late 
times is faster for radiation than for matter dominated branes in agreement with 
6zr ~ and 6zm ~ s~^^^ as found above. 

The interesting case in which the branes possess different equations of state is 
presented in figures |ll] and |12|. Figure [Tl| shows the evolution of the radion when 
the reference brane is radiation dominated and the second brane matter dominated. 



Initially, the trajectories rapidly converge in a similar manner to figures and 10 



however soon the disparity in the expansion rates of the two branes takes effect and 
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Figure 9: The evolution of the interbrane distance z when both branes are radiation 
dominated. The dotted hne is the late time value of the equilibrium position z^- Note the 
way some of the trajectories reach a minimum before approaching their limit: this is due 
to the early time dependence of 

inevitably the interbrane distance reaches zero, corresponding to the branes colliding. 
Figure |1^ shows the opposite case, whereby the reference brane is matter dominated, 
and the second brane is radiation dominated. Again initially the trajectories rapidly 
converge, but since the reference brane is effectively moving faster through the bulk 
than the second brane, the interbrane distance increases indefinitely. 

5.5 Phase Transitions on the Reference Brane 

In this section we investigate the cosmological behaviour of the interbrane distance 
when phase transitions occur on the reference brane. It must be noted that here, 
unlike the previous section we keep the equation of state on the second brane to be 
P2 = —^2 - that of a constant brane tension. 

Here we assume that the reference brane goes through three stages: accelerating, 
radiation dominated and matter dominated. In each stage we ignore all subdomi- 
nant terms in the brane's energy density so that = —1,0,1/3 in each of the 
corresponding phases, leading to, 



A 



A 

\ ao(si)"^ 

ao(si)"^ 
ao(s2)ao(s)^ 



S < Si 
Si < S < S2 
S2 < S 



(5.39) 
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Figure 10: Matter on both branes 




Figure 11: Radiation on the reference brane, matter on the second brane 



where A is a constant and Si and S2 are the times at which the transitions occur. Note 
that here we assume that the transitions are instantaneous. Using this expression 
for rjop and equations ( |5.17| ), (|5.19| ), (|5.20|) and ( ^.21j ) the radion equation ( p.21j ) 
was solved for two cases, each with several initial radion positions distributed evenly 
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Figure 12: Matter on the reference brane, radiation on the second brane 



between z = and z = O.Szhj the results of which are shown in figures |1^ and [14 
Both cases had ri2 = —1.2 and initially r^o = 1.1, and only diff'ered by the choice of 
transition times Si and S2. The first graph, given in figure ^ has Si = 4, S2 — 13 
and shows how during the first phase some of the radion trajectories quickly collide 
with the reference brane, and the rest tend toward Zh as would be expected from the 



analysis of constant tension branes in section |5TT| . When the radiation phase occurs 
however, the surviving solutions quickly move away from the reference brane (in 
these brane based coordinates), before eventually returning to z = 0. The existence 
of the matter phase delays this return further. In the second case given by figure 
the initial phase is twice as long as in the previous situation with si = 8. The 
radiation-matter transition time has been chosen so that matter domination occurs 
at the same energy density as in the first case giving S2 — 17. Doubling Si has an 
interesting effect - the time taken for the inevitable collision between the branes is 
massively delayed, in fact it has been increased roughly by a factor of ten. This 
shows that a relatively short period of stronger infiation on the reference brane than 
on the second brane will delay the return of the second brane indefinitely. The fact 
that in these situations investigated here, the interbrane distance eventually reaches 
zero is of course dependent on the value of ri2 and the limiting value at late time of 
?7o = 1 +?7op. As expected, the behaviour described here and shown in the two figures 
is again in agreement with the conceptual arguments in terms of a bulk observer 
developed in the previous section. 

We now go on to examine the interesting case of two positive tension non-Z2 
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Figure 13: Phase transitions: initially the reference brane has a total constant energy 
density r]o = 1.1, with r]2 = —1.2 on the second brane. At f = 4 and t = 13 the reference 
brane undergoes phase transitions into a radiation dominated universe and then a matter 
dominated universe respectively. 

symmetric branes in a semi-infinite extra dimension. 
5.6 Non-Z2 Symmetry and Weyl Tensor EfFects 

As discussed in section |5]J, it is interesting to examine models where both branes 
possess cosmologically evolving energy densities however, often this requires a second 
brane tension r]2, such that 772 —00 as t 0. This could obviously lead to 
phenomenological difficulties, hence in this section we turn to an interesting model 
in which the reference brane is Z2 symmetric, while the second brane is non-Z2 
symmetric allowing both branes to be of positive tension in a semi-infinite fifth 
dimension. A model of this kind was previously proposed by [58] as an attempt to 
solve the hierarchy problem using only positive tension branes. Several multi-brane 
extensions of this model have been studied [70-72] albeit not for cosmological cases. 

Setting Cq = Ci = k = but assuming C2 7^ 0, means that ho = — 1 and that 
the Friedmann equation for the second brane is given by. 



which can be rewritten using /2 = €2/'^^ as. 



ri2 + ) - 1. (5.41) 
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Figure 14: The initial phase is twice as long as fig [l^ greatly delaying the return of the 
second brane. 



Referring to equations ( |2.1(j| ) and ( |2.2U| ) it can be seen that a solution of the two brane 



system with both branes of positive tension does exist in this situation, provided that, 

— % > V2. (5.42) 

This implies both that /2 < and that the second brane must be dominated by its 
lack of Z2 symmetry, a fact that prevents standard cosmology being realised on this 
brane. We studied the behaviour of the interbrane distance for two different equations 
of state on the second brane. The first was where the second brane possessed a 
constant brane tension r/2 = 1.1 and the results are shown in figure pTS). The reference 
brane had a constant tension of rjQ = 1.01 and the non-Z2 symmetric parameter 
was set to /2 = — 1.0//i^ (in order to keep the radion equation independent of /i) 
and the six initial radion positions were distributed evenly between z = 0.5 and 
z = 2.5. Inspection of equation ( |5.41| ) shows that in this situation the second brane 
will initially expand until a maximum value of 02 (t) is reached which satisfies. 



a2maxiT) = [ 1 . (5.43) 



1/4 

/ir/2(l +r/2). 

Once this is reached /i2 = 0, and subsequently the second brane begins to contract. 
This behaviour helps to explain the trajectories shown in figure |1^: the expanding 
second brane begins to catch up with the reference brane but once 02 — a2max, the 
second brane's motion slows and it quickly falls behind and freezes out at Zh — 2.65. 
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Figure 15: Radion trajectories for tiie two positive tension branes case where = — 1, 
rjQ = 1.01 and r/2 = 1.1 



In the second /2 7^ case investigated, we assumed uj2 = 1/3 and hence that the 
(positive) energy density r]2 was radiation dominated such that 172 = 7r/a^, where 7^ 
is some constant. For each of the six trajectories, we set 7r = 5 and chose a value of 
/2 such that ?72 + ///U?72a^ = —1.1 in order to ensure that the different trajectories 
all start with the same value of = 0.21. The reference brane was assumed to 
have ?7o = 2.0 and the results are shown in figure |TB|. In order to have two positive 
tension branes the inequality ( p.42| ) was maintained and one can therefore see from 
equation (|5.41|) that initially h2 is smaller than hg and the interbrane distance will 
increase, however the r]2 term in equation (|5.41j ) soon becomes negligible and the 
second brane is then entirely dominated by the constant non-Z2 symmetry such that 
/i2 — /2//^^7r ~ 1- Therefore at later times the situation is analogous to that of 
constant tension branes analysed in section |0|: an unstable equilibrium position 
exists as can be seen in figure |1^ and this agrees with the perturbative analysis of 
this situation presented in section |3T2 . 

Obviously there are many other two brane situations that could be investigated 
including both branes having no Z2-symmetry and both branes with non-zero Weyl 
tensor components, however, we leave this to a future investigation. 
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Figure 16: Radion trajectories where the reference brane possesses a constant brane ten- 
sion and the second brane is radiation dominated and non-Z2 symmetric. 

6. Discussion 

In this paper we have greatly extended the work done by Binetruy et al [51] on the 
cosmological radion of the two brane Randall Sundrum model. We started by using 
an elegant method to derive the non-linear equations of motion for the interbrane 
distance 71, and hence found a more general version of the radion equation, while 
sidestepping a lot of the algebra required in [51]. The TZ equation was derived and 
this allowed us to identify the general conditions for an equilibrium radion position. 
The equations of motion were then linearised in order to investigate the nature of the 
equilibrium positions and in the majority of cases these positions were found to be 
unstable, confirming and generalising previous work [47-49]. The radion equations 
were found to be extremely non-linear and it was therefore necessary to numerically 
integrate them in order to examine the radion's behaviour far from the equilibrium 
positions. It should be noted that in some cases no equilibrium exists making a 
numerical analysis even more vital. 

The numerical solutions confirmed that for two constant tension dS branes there 
sometimes existed an unstable equilibrium point, but for two AdS branes that point, 
if it exists would be stable. The cases where the reference brane possessed a more 
realistic equation of state, with standard cosmology taking over at late times were 
also examined, and for a constant second brane tension 772 > 1 it was found that the 
branes inevitably collide. The time taken for this collision to occur was calculated and 
shown to depend upon amoung other things the initial radion position and the nature 
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of the matter on the reference brane, for example, radiation domination implies the 
branes collide sooner than if there was matter domination. It was demonstrated how 
if both branes possess a time dependent energy density then a stable equilibrium 
solution exists provided the branes possess the same equation of state. The effects 
of changes of state were also considered, and it was shown how a relatively short 
burst of inflation could greatly delay the return of the second brane. Two non-Z2 
symmetric scenarios involving two positive tension branes were investigated: one led 
to all trajectories freezing out at Zh, while the other at late times was analogous to 
the the constant tension dS case. It should be noted that these scenarios exist in 
a semi-infinite space y > 0, and assuming /2 < implies that the Schwarzschild 
mass to the right of the second brane C2 < which may lead to problems as the 
second brane will not be shielded from the naked singularity appearing in the bulk. 
This could be avoided by having for example Co — Ci > C2 > i.e. having all the 
Schwarzschild masses greater than zero. 

It was discussed how, by considering the perspective of a bulk observer in a 
static bulk background, one can develop a conceptual understanding of many of 
the features of the radion displayed in this paper. Viewing the behaviour of the 
interbrane distance as a competition between the expansion rates of the branes, 
one can then analyse the branes' Friedmann equations to determine qualitatively 
the radion's behaviour. One may ask whether it would have been better to carry 
out this analysis in the bulk based coordinate system as opposed to a brane based 
one. However, we intend to extend this radion analysis to systems containing more 
complicated bulk matter than just a cosmological constant, namely scalar and dilaton 
fields. It is then not possible to use a time independent bulk based coordinate system 
since Birkhoff's theorem (the generalisation of Gauss' theorem) can no longer be 
applied. Given the elegance and ease of deriving the radion equation in the brane 
based system, the need to follow and extend the work done by [51], and the inevitable 
numerical calculation whatever coordinate system was used, it was felt that the brane 
based system was the most appropriate for our purposes. 

Finally, it should be mentioned how special a case the cosmological radion is. 
Due to the amount of symmetry assumed, the trajectories of the branes are entirely 
determined by the knowledge of their Friedmann equations and hence their expansion 
rate, a fact which allowed us to efficiently derive the radion equation, and was not 
realised by [51] who had to use a more cumbersome method. In essence the branes 
do not 'feel' each other, a fact that is altered if one considers other fields in the bulk, 
or general perturbations to the cosmological background. 

Acknowledgements: We would like to thank Pierre Binetruy and Carsten van de 
Bruck for useful discussions. This work is supported in part by PPARC. 



36 



References 



V. A. Rubakov and M. E. Shaposnikov, Do we live inside a domain wall?, Phys. Lett. 
B125 (1983) 136-138. 

K. Akama, An Early Proposal of "Brane World", Led. Notes. Phys 176 (1982) 



267-271, [lhep-th/0001113 |. 



I. Antoniadis, A possible new dimension at a few TeV, Phys. Lett. B246 (1990) 377. 

M. Visser, An exotic class of Kaluza-Klein models, Phys. Lett. B159 (1985) 22, 
[|hep-th/99 100931] . 



N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, The Hierarchy Problem and New 
Dimensions at a Millimeter, Phys. Lett. B429 (1998) 263-272, ||hep-ph/9803315| . 



N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, Phenomenology, Astrophysics and 
Cosmology of Theories with Sub- Millimeter Dimensions and Te V Scale Quantum 
Gravity, Phys. Rev. D59 (1999) ||hep-ph/9807344|| . 



I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, New Dimensions at a 
Millimeter to a Fermi and Superstrings at a TeV, Phys. Lett. B436 (1998) 
[|hep-ph/9804398|] . 



R. Sundrum, Effective Field Theory for a Three-Brane Universe, Phys. Rev. D59 
(1999) 85009, [hep-ph/980547lj] . 

L. Randall and R. Sundrum, An Alternative to Compactification, Phys. Rev. Lett. 83 



(1999) 4690-4693, ||hep-th/9906064|] 



L. Randall and R. Sundrum, A Large Mass Hierarchy from a Small Extra Dimension, 



Phys. Rev. Lett. 83 (1999) 3370-3373, |tiep-ph/9905221|. 



P. Binetruy, C. Deffayet, and D. Langlois, Non-conventional cosmology from a 



brane-universe, Nucl. Phys. B565 (2000) 269-287, ||hep-th/9905012|| 



P. Binetruy, C. Deffayet, U. Ellwanger, and D. Langlois, Brane cosmological evolution 
in a bulk with cosmological constant, Phys. Lett. (2000) 285-291, | hep-th/991021£ l. 



C. Csaki, M. Graesser, C. Kolda, and J. Terning, Cosmology of One Extra Dimension 
with Localized Gravity, Phys. Lett. B462 (1999) 34-40, [|hep-ph/990"6513| . 



J. M. Cline, C. Grojean, and G. Servant, Cosmological Expansion in the Presence of 
an Extra Dimension, Phys. Rev. Lett. 83 (1999) 4245, [|hep-ph/9906523|| . 

P. Kraus, Dynamics of Anti-de Sitter Domain Walls, JHEP 9912 (1999) Oil, 
||hep-th/9910149|] . 

S. Mukohyama, Gauge-invariant gravitational perturbations of maximally symmetric 



spacetimes, Phys. Rev. D62 (2000) 084015, [|hep-th/0004067 | 



37 



H. Kodama, A. Ishibashi, and O. Seto, Brane World Cosmology - Gauge-Invariant 
Formalism for Perturbation, Phys. Rev. D62 (2000) 064022, [|hep-th/0004160 |. 



D. Langlois, Brane cosmological perturbations, Phys. Rev. D62 (2000) 126012, 



||hep-th/0005025|| . 



C. van de Bruck, M. Dorca, R. Brandenberger, and A. Lukas, Cosmological 
Perturbations in Brane-World Theories: Formalism, Phys. Rev. D62 (2000) 123515, 



||hep-th/0005032|l . 

R. Maartens, D. Wands, B. A. Bassett, and I. P. C. Heard, Chaotic inflation on the 



brane, Phys. Rev. D62 (2000) 041301, [|hep-ph/9912464 |. 

D. Langlois, R. Maartens, and D. Wands, Gravitational waves from inflation on the 



brane, Phys. Lett. B489 (2000) 259-267, [lhep-th/0006007|| . 



S. C. Davis, W. B. Perkins, A. C. Davis, and I. Vernon, Cosmological Phase 
Transitions in a Brane World, Phys. Rev. D63 (2001) 083518, | |hep-ph/00l"222^ ]. 



A. Chamblin and G. W. Gibbons, Supergravity on the Brane, Phys. Rev. Lett. 84 
(2000) 1090-1093, ||hep-th/9909130(1 . 



T. Gherghetta and A. Pomarol, Bulk Fields and Supersymmetry in a Slice of AdS, 
Nucl. Phys. B586 (2000) 141-162, [|hep-ph/0003T29| . 



R. Altendorfer, J. Jonathan Bagger, and D. Nemeschansky, Super symmetric 



Randall- Sundrum Scenario, Phys. Rev. D63 (2001) 125025, |hep-th/0003117|l . 
M. Duff, J. T. Liu, and K. Stelle, A super symmetric Type IIB Randall- Sundrum 



realization, J. Math. Phys. 42 (2001) 3027-3047, |hep-th/000712C|] . 

E. Bergshoeff, R. Kallosh, and A. Van Proeyen, Supersymmetry in Singular Spaces, 



JHEP 0010 (2000) 033, ||hep-th/0007044|] . 

C. S. Chan, P. L. Paul, and H. Verlinde, A Note on Warped String Compactification, 



Nucl. Phys. B581 (2000) 156-164, |lhep-th/0003236|| . 

M. Cvetic, M. Duff, J. T. Liu, H. Lu, C. Pope, and K. Stelle, Randall- Sundrum 



Brane Tensions, Nucl./ hys. B605 (2001) 141-158, [|hep-th/00lTl67 |. 

P. Brax and A. C. Davis, Cosmological Solutions of Supergravity in Singular Spaces, 



Phys. Lett. B497 (2001) 289-295, [|hep-th/0011045 |. 

P. Brax and A. C. Davis, Cosmological Evolution on Self-Tuned Branes and the 



Cosmological Constant, JHEP 0105 (2001) 007, ||hep-th/0104023|l . 

P. Brax and A. C. Davis, On Brane Cosmology and Naked Singularities, Phys. Lett. 



B513 (2001) 156-162, [ |hep-th/0105269 | 



38 



p. Brax, C. van de Bruck, and A. C. Davis, Bmne-World Cosmology, Bulk Scalars 



and Perturbations, JHEP 0110 (2001) 026, [|hep-th/0108215 | 



P. Brax, C. van de Bruck, A. C. Davis, and C. Rhodes, Wave function of the radion 
with a bulk scalar field, Phys. Lett. B531 (2002) 135-142, ||hep-th/020119l|| . 



P. Brax, C. van de Bruck, A. C. Davis, and C. Rhodes, Four-dimensional Gravity 



from Singular Spaces, Phys. Rev. D65 (2002) 121501, ||hep-th/020213^ . 

P. Brax, C. van de Bruck, A. C. Davis, and C. Rhodes, Cosmological Evolution of 



Brane World Moduli, Phys. Rev. D67 (2003) 023512, [|hep-th/0209158 | 



R. Gregory, Nonsingular global string compactifications, Phys. Rev. Lett. 84 (2000) 
2564-2567, ||hep-th/9911015|| . 

T. Gherghetta and M. Shaposhnikov, Localizing Gravity on a String-Like Defect in 
Six Dimensions, Phys. Rev. Lett. 85 (2000) 240-243, ||hep-th/0004014|| . 



I. Olasagasti and A. Vilenkin, Gravity of higher-dimensional global defects, Phys. 



Rev. D62 (2000) 044014, ||hep-th/0003300 |. 



T. Gherghetta, E. Roessl, and M. Shaposhnikov, Living Inside a Hedgehog: 
Higher- dimensional Solutions that Localize Gravity, Phys. Lett. B491 (2000) 
353-361, [|hep-th/000"625l |. 



P. Bostock, R. Gregory, I. Navarro, and J. Santiago, Einstein gravity on the 



codimension 2 brane?, |hep-th/0311074 



W. D. Goldberger and M. B. Wise, Modulus Stabilization with Bulk Fields, Phys. 



Rev. Lett. 83 (1999) 4922-4925, ||hep-ph/9907447|| . 



C. Csaki, M. Graesser, L. Randall, and J. Terning, Cosmology of Brane Models with 
Radion Stabilization, Phys. Rev. D62 (2000) 045015, |hep-ph/9911406|] . 



W. D. Goldberger and M. B. Wise, Phenomenology of a Stabilized Modulus, Phys. 



Lett. B475 (2000) 275-279, |lhep-ph/9911457|| . 

T. Tanaka and X. Montes, Gravity in the brane-world for two-branes model with 



stabilized modulus, Nucl. Phys. B582 (2000) 259-276, | tiep-th/0001092 1 . 



S. Mukohyama and L. Kofman, Brane Gravity at Low Energy, Phys. Rev. D65 
(2002) 124025, ||hep-th/0112115|| . 



C. Charmousis, R. Gregory, and V. A. Rubakov, Wave function of the radion in a 
brane world, Phys. Rev. D62 (2000) 067505, ||hep-th/9912160|| . 



U. Gen and M. Sasaki, Radion on the de Sitter brane. Prog. Theor. Phys. 105 (2001) 
591-606, fer-qc/001107^ . 



39 



A. Karch and L. Randall, Locally Localized Gravity, JHEP 0105 (2001) 008, 
||hep-th/0011156|] . 



C. Csaki, M. Graesser, L. Randall, and J. Terning, Cosmology of Brane Models with 
Radion Stabilization, Phys. Rev. D62 (2000) 045015, |hep-ph/99l"l40^ . 

P. Binetruy, C. Deffayet, and D. Langlois, The Radion in Brane Cosmology, Nucl. 



Phys. B615 (2001) 219-236, |hep-th/0101234|. 



A.-C. Davis, C. Rhodes, and I. Vernon, Branes on the Horizon, JHEP 0111 (2001) 



015, |lhep-ph/0107250|| 



A.-C. Davis, I. Vernon, S. C. Davis, and W. B. Perkins, Brane World Cosmology 
Without the Z2 Symmetry, Phys. Lett. B504 (2000) 254-261, ||hep-th/0008132| . 



W. Israel, Singular Hypersurfaces and thin shells in general relativity, Nuovo. Cim. 
B44S10 (1966) 1. 

G. Darmois, Memorial des Sciences Mathematique XXV, Gauthier - Villars, Paris 
(1927). 

C. Misner, K. Thorne, and J. Wheeler, Gravitation, Freeman (1982). 

R. Wald, General Relativity, University of Chicago Press (1988). 

J. Lykken and L. Randall, The Shape of Gravity, JHEP 0006 (2000) 014, 



||hep-th/9908076|| . 



S. Mukohyama, T. Shiromizu, and K.-i. Maeda, Global structure of exact cosmological 



solutions in the brane world, Phys. Rev. D62 (2000) 024028, ||hep-th/9912287|| . 



N. Kaloper, Bent Domain Walls as Braneworlds, Phys. Rev. D60 (1999) 123506, 



iep-th/9905210|l . 



C. Csaki, J. Erlich, and C. Grojean, Gravitational Lorentz Violations and Adjustment 
of the Cosmological Constant in Asymmetrically Warped Spacetimes, Nucl. Phys. 
B604 (2001) 312-342, [|iep-th/00l"2l43 |. 



J. Khoury, B. A. Ovrut, P. J. Steinhardt, and N. Turok, The Ekpyrotic Universe: 
Colliding Branes and the Origin of the Hot Big Bang, Phys. Rev. D64 (2001) 123522, 
||hep-th/0103239|| . 



R. Kallosh, L. Kofman, and A. Linde, Pyrotechnic Universe, Phys. Rev. D64 (2001) 



123523, ||hep-th/0104073|| 



J. Khoury, B. A. Ovrut, P. J. Steinhardt, and N. Turok, Density Perturbations in the 
Ekpyrotic Scenario, Phys. Rev. D66 (2002) 046005, ||hep-th/0109050|| . 



P. J. Steinhardt and N. Turok, Is Vacuum Decay Significant in Ekpyrotic and Cyclic 



Models?, astro-ph/0112537 



40 



[66] A. J. Tolley and N. Turok, Quantum Fields in a Big Crunch/Big Bang Spacetime, 
Phys. Rev. D66 (2002) 106005, [tiep-th/020409l |. 

[67] A. J. Tolley and N. Turok, Cosmological Perturbations in a Big Crunch/Big Bang 
Space-time, |hep-th/0306109| . 

[68] P. J. Steinhardt and N. Turok, Cosmic Evolution in a Cyclic Universe, Phys. Rev. 
D65 (2002) 126003, [ |liep-th/01lT098 |. 



[69] P. J. Steinhardt and N. Turok, A Cyclic Model of the Universe, |hep-th/0111030 . 

[70] S. Mouslopoulos and A. Papazoglou, "+-+" Brane Model Phenomenology, JHEP 
0011 (2000) 018, |hep-ph/0003207[| . 

[71] I. I. Kogan, S. Mouslopoulos, A. Papazoglou, G. G. Ross, and J. Santiago, A Three 
three-hrane Universe: New Phenomenology for the New Millennium?, Nucl. Phys. 
B584 (2000) 313-328, [>iep-ph/9912552 |. 

[72] I. I. Kogan and G. G. Ross, Brane Universe and Multigravity: Modification of gravity 
at large and small distances, Phys. Lett. B485 (2000) 255-262, [ fiep-th/0003074 |. 



41 



